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Abstract 

-I— > 

We study the propagation of acoustic waves in a fluid that is contained in a 
, thin two-dimensional tube, and that it is moving with a velocity profile that only 

depends on the transversal coordinate of the tube. The governing equations are the 
Galbrun equations, or, equivalently, the linearized Euler equations. We analyze 
the approximate model that was recently derived by Bonnet-Bendhia, Durufle 
^ and Joly to describe the propagation of the acoustic waves in the limit when the 

width of the tube goes to zero. We study this model for strictly monotonic stable 
|>! velocity profiles. We prove that the equations of the model of Bonnet-Bendhia, 

Durufle and Joly are well posed, i.e., that there is a unique global solution, and 
that the solution depends continuously on the initial data. Moreover, we prove 
that for smooth profiles the solution grows at most as t 3 as t — > oo, and that for 
piecewise linear profiles it grows at most as t 4 . This establishes the stability of the 
model in a weak sense. These results are obtained constructing a quasi-explicit 
representation of the solution. Our quasi-explicit representation gives a physical 
interpretation of the propagation of acoustic waves in the fluid and it provides an 
efficient way to compute numerically the solution. 



*AMS 2000 classification: 35Q72; 35Q35; 35B35; 35C15. 

^On leave of absence from Departamento de Metodos Matematicos y Numericos. Instituto de Investi- 
gaciones en Matematicas Aplicadas y en Sistemas. Universidad Nacional Autonoma de Mexico. Apartado 
Postal 20-726, Mexico DF 01000. Fellow of the Sistema Nacional de Investigadores. 



1 



1 Introduction 



In this paper we study the following initial value problem, 



' find u(x,y,t) : R x [-1, l]xl + ^ R, 
^dt +M{y) di )2u ~2d^ udy = 0, (x,y)eRx[-l,l), t > 0, 

(1.1) 

u(x,y,0) = u°(x,y), (x,y) eRx [-1,1], 

— (x,y,0) = u l (x,y), (x,y) ERx [-1,1], 

where M G L°°([— 1, 1]) is a real-valued function, and the initial data u° G Ly(H™ +1 ) and 
u 1 E L 2 y (H™),m = 0, 1, ■ ■ ■ are given, with 

Ll(H™):=L 2 ([-l,l);H™(R)). (1.2) 

We look for solutions in "natural energy spaces" of the form 

u G C° (R + ; Ly(H™)) n C 1 (R+; L^H™' 1 )) , m = 0, 1, • • • . (1.3) 

This mathematical model has been obtained by Bonnet-Bendhia, Durufle and Joly |3] 
as an approximation in the study of a problem of aero-acoustics. They considered the 
propagation of acoustic waves in two dimensions in a fluid that is contained in a thin 
tube. The function M{y) describes, in normalized coordinates, the lateral variations of 
the velocity of the fluid. Actually, x is the coordinate along the axis of the tube and y 
is the transversal coordinate. The velocity of the fluid is directed along x and it is given 



by M(y) at the point (x,y) of the tube. The model (1.1) was obtained in [3] by a formal 
asymptotic expansion on the width of the tube of the solution to Galbrun equations [5], 
that are equivalent to the linearized Euler equations. Physically, the validity of our model 
requires that the transverse dimension of the tube is small with respect to the wavelength 
but no too small to justify the fact of neglecting the viscosity effects : in this sense, this 



model can be seen as a low frequency model. Another potential application of this study is 
the construction of effective boundary conditions - also called lining models - to take into 
account the modeling of boundary layers in aeroacoustics, which is a quite delicate issue 
from mathematical and numerical points of view (see for instance [3]). 



In |3j the analysis of equation (1.1) was reduced to a one dimensional problem by 



Fourier transform along x, taking advantage of the fact that the velocity profile is only a 



function of the transversal coordinate y. This reduces the study of the solutions to (1.1) 
to the spectral theory of a non-local operator, A, that acts only on the transversal variable 
ye [-1,1]: Ae C(L 2 [-1, l]) 2 ). More precisely, setting 

U(x,y,t) := (u{x,y,t),(d t + M(y)d x )u{x,y,t)Y 

and denoting by U(k,y,t) the Fourier transform along x of U(x,y,t) one has formally 

^-(k,;t)+ikAU(k,-,t)=0 U(k, -,t) = e~ lkAt U(k, -,0). (1.4) 

It was shown in [3] that on spite of its apparent simplicity this problem has rather surprising 
properties. As the operator A is not normal, it can have complex (non-real) eigenvalues. 



It was pointed out in [3J that a necessary condition in order that the problem (1.1) is well 
posed and stable - in the sense that the solutions do not grow exponentially in time- is that 
all the eigenvalues of A are real. These issues were studied in detail in [3J. General prop- 
erties of A were obtained and the general structure of the point and continuous spectrum 
were analyzed. Moreover, several results on the existence and on the absence of complex 
eigenvalues of A were given. These results were illustrated numerically. Furthermore, it 
was conjectured in [3] that the condition that A has no complex eigenvalues is also suffi- 



cient for the stability of the problem (1.1). Unfortunately, it does not seem that such a 
result can easily be deduced from the standard semi-group theory. 

That is why, in this paper, we take a point of view that is slightly different from the one 



of [3j. Actually, we directly solve equations (1.1) by Fourier-Laplace transform. As usual 



to prove that the solution does not grow exponentially in time it is necessary to deform the 
integration contour of the inverse Laplace transform to the real axis. This requires that 



the norming coefficient, N(X), in the inverse Laplace transform (see (3.7)) has no complex 
poles. It turns out that the poles of N(X) in C \ Range M are precisely the eigenvalues of 
the operator A of [3j in C \ Range M, what means that we obtain, as to be expected, the 
same stability condition as in [3J. Moreover, N(X) has, in general, a cut on Range M and, 
furthermore, it can blow up as we approach the cut from above and from below, what makes 
the issue of deforming the contour to the real axis on both sides of the cut quite delicate. 
Although we think that our approach is quite general, it is difficult to state general results. 
The method that we present here to construct a quasi-explicit representation of the solution 
and to prove stability can be applied to general profiles, provided that one can analyze the 
limiting values of the norming factor N(X) as A tends to [M_,M + ] from above and from 
below. In this paper, we develop the well-posedness theory in two cases. First, we consider 
strictly monotonic and convex (or concave), smooth profiles. In this case, we prove that 
the norming factors N(X) of our profiles have continuous limiting values as we approach 
the cut from above and from below, that are different. We use this result to prove that the 



problem (1.1) has a unique solution in the natural class of fields (1.3). Moreover, we prove 
that the solutions grow at most as t 3 as t — > oo. Then, we consider the case of piecewise 
linear profiles for which the structure of N(X) is quite different. In our mind, the interest 
of this class of profiles (already considered in [3] as a theoretical tool) is to provide a safe 
way to get numerical approximations of the solution. We also prove that for these profiles 



the problem (1.1) has a unique solution in the natural class of fields (1.3). However, in this 
case we prove that the solutions grow at most as t 4 as t — > oo. These results establish that 
the model is stable in a weak sense, and prove the conjecture of [3 J for our class of profiles. 
For both cases of profiles our method, that is based in the Fourier-Laplace transform, leads 
to a quasi-explicit representation of the solution, that gives a physical interpretation of 
the propagation of acoustic waves in the fluid, and it provides an efficient way to compute 
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numerically the solution. The numerical results are presented in |6]. 

The paper is organized as follow. In Section 2 we briefly state, for the reader's conve- 



nience, a derivation of the approximate model (1.1) from the Galbrun equations [5] that 
is slightly different from the one in [3]. In Section 3 we develop the well-posedness theory 
and we prove our results on the existence and uniqueness of solutions in the space of fields 
(T3J) and on the continuous dependence on the initial data for the two classes of profiles 



(Subsections 3.3.1 and 3.3.2, respectively). Finally, in Section 4, we obtain a quasi-explicit 



representation of the solutions to (1.1), of which we give a physical interpretation. 



2 Derivation of the quasi- ID model 



We model this flow by means of the equations of Galbrun [5] where the unknown (u e , v e ), 
that are functions of (x,y,t), are the components of the Lagrange displacement of the 
fluid. The governing equations are, 



( ( d lf , , ()\ 2 d ( du £ d\ £ \ „ , . _ 



(2.1) 



d ( du £ d\ £ 

+ 



0, (x,y)eQ £ , t>0, 



dy \ dx dy 

with the boundary condition that on the walls of the tube, y = ± e, the normal component 
of the displacement is zero, 



v e (x, ± e, t) = 0, xeR, t>0. 



(2.2) 



To determine formally the approximate model (see [3j) let us first perform a change of scale 
in y (y := y / s) to work in a fixed domain and use a scaling of the unknown v £ (which we 
expect to tend to since this function is defined on a thin domain and is imposed to on 
the boundary of this thin domain). Let us denote, 

u £ (x,y,t) := u £ (x,ey,t), (x, y) :=£ R X (—1, 1), t > 0, 
v £ (x,y,t) := v e (x,ey,t)/e, (x, y) <E R X (-1, 1), t > 0. 



(2.3) 



The equations in (u e ,v £ ) are given by. 
+ M(y) -)u 



dt 



\ 2 d ( du £ dv £ 



dx \ dx dy 



= 0, (x,y) 6 R x (-1,1), t>0, 



with the boundary conditions, 

v e (x,±l,t) = 0, iGl, t>0. 



(2.4) 



(2.5) 



We obtain the limit model by dropping the e 2 term in the second equation of (2.4) : the 

mi 
( <d 



formal limit (u,v) of (u £ ,v £ ) satisfies 



Q \ 2 d ( du dv 
dt + M ' y ' dx) u ~8x \dz + dy) = °' (*•») ^ » x (-1,1), t > 0, 



d f du dv s ^ 



(ar,j/) Glx (-1,1), *>0. 



Then, there is a function of x,t only, that we denote by p(x,t), such that, 

/du dv s 
\dx dy, 

Integrating in y and as v(x, ± l,t) = 0, we prove that, 

1 Z" 1 9m 



p(x,t) 



2 y_! &r 



[x,y,t) dy. 



The first equation (2.4) writes 



d „ w n <9 \ 2 d f du dv 



dt 



dx \dx dy 



d_ 

dt 



+ M(y) -)« + 



<9 ^ 2 9p 
<9x 



0. 



Finally, replacing p by its value given by (2.7) we obtain equation (1.1). 



(2.6) 



(2.7) 



3 Well-posedness of the quasi- ID model 



3.1 Preliminary material 



In what follows we assume that M(y) is either increasing or decreasing and we denote, 
M+ := sup M(y), M_ := inf M(y). 

ye (-1,1) ye (-1,1) 
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We write (ll.ll) as follows 



(| + M(y) |) 2 «(x,,,t) 



<9x 2 



[a(u)](z,*), 



where v — > a(u) is the averaging operator in the y— direction: 



1 
2 



(3.1) 



(3.2) 



From (3.2), it is clear that, if a(u) is known, (3.1) is a simple transport square equation 
along the axis of the tube for each fixed y, with a y-dependent transport velocity that is 
given by the profile M(y). We solve this transport equation explicitly. In this way we 



obtain a quasi-explicit representation of the unique solution to (1.1). Therefore, we would 



like to get an equation for the average value a(w) alone. As we shall see, such an equation 
in a(u) is obtained by applying the Fourier transform in x and the Laplace transform in 



t. Then, the well-posedness of (1.1) is reduced to showing that it is possible to apply the 



inverse Fourier-Laplace transform, and also to obtain a priori estimates for a(w). 

As usual, we define the Fourier transform as an unitary operator on L 2 (IR), defined for 
/ G L\R) as 



(.F/)(fc) = /(*) := 



1 



-ikx 



f(x) dx. 



After Fourier transform equation (3.1) becomes (this is nothing but (1.4) written is scalar 
second order form), 







(— + M(y)ik) u = -k 2 a(u)(k,t), keR, t > 0. 



(3.3) 



We take a point of view that is slightly different from the one of [3j. Instead of writing 



equation (3.3) as an evolution problem for a first order system, and using semigroup theory 



to reduce the problem to the spectral analysis of the generator A(k) = k A of the system, 



as it was done in [3], we directly solve equation (3.3) by Laplace transform in time. 



We denote by u the Laplace transform in time of u, 



u(k,y,u):= / e luJt u(k,y,t)dt, Slu > 0. 
'o 



(3.4) 



REMARK 3.1. One easily checks that, for each k e E, ( |3.3P admits a unique solution 
u(k,-,t) which increases in time at most exponentially as (because A is bounded see 

expS|/c|t with B := \\A\\ > 0, 
so that u(k, -,t) has a well defined Laplace transform as soon as > B\k\. 



After Laplace transform, (3.3) becomes 



(oj - M(y) k) 2 u + k 2 Si(u)(k, t)=u l - i(u - 2kM(y))u°, (3.5) 



where u(x,y, 0) = u°(x,y), and —u(x,y,0) = Dividing both sides of (3.5) by 



(u — kM(y)) 2 and taking the average over y of both sides we prove that, 

a(w) = -2N a ( [u 1 - i(u - 2kM{y))u°] [u - kM(y)] ~ 2 ) , (3.6) 
where N is the norming factor, 

iV(A):=(2-F(A))- 1 , (3.7) 

with, 

F(X):= [ (M(y)-\y 2 dy, XeC\[M.,M+]. (3.8) 



i 



Inverting the Laplace transform and changing the integration variable u to A = \r + i\j := 
u)/k, we obtain the following representation for a(u), 

&(u)(k,t) = ao(k,t) + SLi(k,t) (3.9) 

where &o(k,t) and ki(k,t) are, respective, the contributions of the Cauchy data u° and u 1 
and are given by: 

7 — j P 

a Lj (k,t) = -— ^ e- ikXt N(X) ai(f j (;X)u j (k,-))dX, (3.10) 

where A/ > for k > and A/ < for k < with |A/| large enough. Furthermore, 

f (y, A) := iM{y) (A - M(y)y 2 - i (A - M(y))-\ A(y, A) := (A - M(y))- 2 . (3.11) 



For proving the well-posedness of (1.1), our objective is to derive a priori estimates for 



ao(k,t) and a.i(k,t) as functions of k and t. Such estimates can not be deduced directly 



from ( 3.10, 3.11) as if A is not real, the exponent in the right-hand side of (3.10) will blow 
up as k — > ±oo, unless the initial data decays very fast as \k\ — > oo, and we will not be able 
to invert the Fourier transform to compute [a(it)] (x,t). It is now imperative to transform 



the expressions (3.10, 3.11) using complex integration techniques (contour deformation). 
This is the object of the next subsection. 



3.2 A new expression for ao(k,t) and &i(k,t) 

We consider the case k > 0. For k < the results are obtained in the same way, with 
obvious changes. 



According to (3.10) and (3.11), we simply need to compute the integrals 

N(X) 



h(kt,y) := I < 



-ikXt 



(A - M(y)) 



e+i 



dX, A/ > 0, 



0,1. 



(3.12) 



Indeed, we immediately obtain from ( 3.10, 3.11) 



ao(M) = — a([M(.)Ji(fa,0 -I (kt,-)] u°(k,-) 
*i(M) = ^ a^A^Ou 1 ^,-))- 



(3.13) 



(3.14) 



We are going to compute the integrals (3.12) by using complex integration methods. For 



this, we need to describe in detail the analyticity properties of N(X). From now on, we 
shall assume that, 

(Hi) M is a continuous function of y. 

We set 



M_ = min M(y), M 4 



max M(y). 



(3.15) 
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From formula (3.8), it is clear that F is analytic in C\ [M_, M+] so that N is meromorphic 
in C \ [M_, M+]. Observe that as 

lim N(X) = \ 



the poles of N(X) are contained in a bounded set of C. Since N(X) = N(X), we deduce 
that these poles are equally distributed with respect to the real axis. At this point, we 
need to make a fundamental assumption 

(Hs) The function N(X) has no complex (i.e. non real) poles. 

This assumption is clearly related to the velocity profile M. The paper |3] gives explicit 
conditions on M that ensure that (Hs) is satisfied. We shall give another proof of such a 



result in Lemma 3.3 The paper [3] also gives examples of profiles for which (Hs) does not 
hold. We shall say that a profile M is stable when the assumption (Hs) is satisfied. This 
denomination is justified by the fact that one easily proves that if (Hs) does not hold, the 
Cauchy problem is strongly ill-posed (see [3] - the complex poles of N(X) are precisely the 
complex eigenvalues of the operator A of [3]). The assumption (Hs) is thus a necessary 



condition for the well-posedness of (1.1). Our conjecture is that this condition is also 
sufficient. The object of this paper is to show that this conjecture is true under additional 
technical assumptions on M. 

Concerning the real poles of N(X), we note (see again [3] for details) that N(X) has at most 
one pole in (— oo, M_) that we denote by A_, if it exists, and at most one pole in (M+, oo) 
that we denote by A + if it exists, and that these poles are simple. This is immediate since, 

d f 1 1 

F(X) = -2 / — — — dX 



OX " ' J_ x (M(y) - A) 3 

is positive for A > M + and is negative for A < M_. Furthermore, the residues of N(X) at 
the poles are given by, 

Res N(X ± ) = j=±- y (3.16) 
10 



We make the additional assumption, 

('Hrl) M(y) takes the value M± at a point where it has a right or a left derivative. 

One easily sees (see also [3]) that in this case F(X) blows up to +00 when A approaches 
M±, from real values, which ensures the existence of A±. 

We now present our computations of li in the case k > 0. As mentioned above, for 
k < the results are obtained in the same way. First, thanks to (Hs), we can deform 
the line Im A = A/ into the contour that coincides with the real axis outside a 5- 
neighborhood of the poles A± and the cut [M_,M + ] (see Figure [T]) to obtain (note that 
N(X)/ (A — M(y)Y +1 decays to when \lZe A| goes to infinity): 



h{kt,y) = I e~ ikXt - N [ X \, +1 d\. (3.17) 



(A - M(y)y 



Im X 



Tm A = Aj 



A_ 



A + Ue X 



C + 

^6 



Figure 1: The integration contour C s 



Let C s be the contour deduced from Cf by symmetry with respect to the real axis (see 
Figure g). 

We claim that, 

' - NW dA = 0. (3.18) 



-ikXt 



(A - M(y)) 
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Im A 



fte A 



Figure 2: The integration contour C$ 



This is obtained by considering the closed contour C R s of Figure |3| inside which the inte- 
grand is analytic. 



Im A 




Figure 3: The integration contour Cr,, 



Tie A 



Thus, by Cauchy's theorem 



-ikXt 



N(X) 



(A - M(y)) 



t+i 



d\ = 0. 



(3.19) 



Then, (3.18) is obtained from (3.19) by passing to the limit when R — > +oo: the contri- 



bution of the integral along the semi-circle of radius R vanishes because it is inside the 
good complex half-space when k > (for the exponential term) and because the function 



iV(A)/(A - M(y)Y +1 decays at infinity 



Finally adding (3.17) and (3.18), we obtain 



ikXt 



N(X) 



h(kt,y) = I e 

Ics (A - M(y)) H 



j d\, 



(3.20) 
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where Cs is the closed contour of Figure [4} 



Tm X 




Tie A 



Figure 4: The integration contours Cs and Ss 



By the theorem of residues we also have: 



I t (kt, y) = Ii, p {kt, y) + h )C (kt, y), 



1,2, 



(3.21) 



where 



h, p {kt,y) : = -2m ^ 



-ik\±t 



ResiV(A d 



L c (kt,y) :-- 



± 

-ikXt 



N(X) 



(3.22) 



is* (A - M(y)) l+1 

where Ss is a closed curve that collapses to the segment [M_,M + ] when 5 goes to (see 

Figure lil again). 



To get a more exploitable formula for Ii )C (kt,y) , we would like to take the limit of 
the expression in the right-hand side when 5 tends to 0. For this we need some additional 
information on the behaviour of A^(A) when A approaches the segment [M_,M + ]. In the 
following we consider two different situations leading to quite different behaviours of N(X) 
when A approaches [M_ , M + ] . 

• The case of a class of smooth profiles with constant monotonicity and convexity (see 



Subsection 3.3.1). This is a physically relevant case for which N(X) has two distinct 
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limits on [M_,M + ] depending on whether one approaches the segment from above 
or from below. In this case, the function N(X) presents a cut along [M_, M + ). 

• The case (already considered in [3]) of a class of piecewise linear profiles (see Sub- 
section 3.3.2), typically with fixed convexity . This case is interesting to provide a 
safe procedure for numerical approximations. In this situation, the function N(X) is 
meromorphic (it is a rational fraction) with real poles distributed in [M_,M + ], as 
well as poles at A±. 

3.3 Well-posedness results 
3.3.1 For a class of smooth profiles 

We denote by C^[-l,l],n = 0, 1, ■ ■ ■ ,0 < 7 < 1, the space of real-valued functions 
defined on [—1, 1] that are n— times differentiable and with the derivative of order n Holder 
continuous with exponent 7. 

DEFINITION 3.2. If I denotes an integer > 2 and 7 a real number in ]0, 1[ . we shall 
say that a profile M(y) belongs to the class C 1 ^ 1 if it fulfills the following conditions: 

Me^[-l,l], M\y)?0, M"(y)^0, Vye[-l,l]. (3.23) 

In other words, a profile in the class C 1 ^ is smooth, strictly monotonous and has a fixed 
(strict) convexity. 

LEMMA 3.3. Suppose that M(y) belongs to the class C l f for some I > 2 and 7 e]0, 1[. 
Then, it is stable in the sense that it satisfies the condition (Hs)- 

Proof: Without loss of generality, we can assume that M is increasing (changing M into 
— M simply exchanges the signs of the roots of F(X) = 2) and that M_( = M{— 1)) = 
(translating M simply translates the solutions of F(X) = 2 parallel to the real axis). Since, 
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by continuity, M'(y) is bounded from below by a strictly positive constant, we can write 

1 M'(y) dy f 1 d f 1 \ dy 



F(X) 1 i (M(y) - XT M'(y) 
After integration by parts, we get 

1 M"(y) 



i - 



x dy \(M(y)-X)j M'(yY 



F(X) 7-i (M(y) - A) M'{yY 
That is to say, for A = A_r + iA/, 

- 1 M(y) - X R M"(y) 



dy 



1 



Tie F(X) 



_x |M(y)-A| 2 M'(y) 2 ^ + 



Im F(A) = -A/ 



.! \M{y)-X\ 2 M'(y) 2 
If -F(A) = 2 with A/ 7^ 0, we have in particular 

1 M"{y) 



dy + 



(M(y) - A) M'(y) 

M(y) - X R 
l\M(y)-X)\ 2 M>(y) 

1 



1 

-1 



1 

-1 



|M(y)-A)| 2 M'(y) 



-1 



dy + 



1 



0. 



-1 



(3.24) 



U \M (y) - X\ 2 M'(y) 2 » V\M{y)-X)\ 2 M'{y). 
Substituting this in the right hand side of TZe F(X) (note that the above term is the 
multiplicative factor of A#), we get (since M(— 1) = and setting M' ± = M'(±l) ) 



Tie F(X) 



M{y) M"{y) 



dy - 



Mi 



Li \M(y) - A| 2 M'(y) 2 a \M+ - A| 2 M\ 
Next, we claim that TZe F(X) < so that , F(X) = 2 is impossible. We distinguish two 
cases: 



(i) M"(y) > 0. In this case, TZe F(X) is clearly the sum of two negative terms. 



(ii) M"(y) < 0. In this case, we can write 
•' M{y) M"{y) 



_! |M(y)-A| 2 M'(y) 2 



dy < 



Mi 



M"{y) 
_! \M{y) - A| 2 M'{y) 2 



|M(y)-A)| 2 M'(y) 



The second inequality being deduced from (3.24). Then, 



Tie F(X) < 



M 4 



Mi 



\M(y) - X)\ 2 M'(y)\-i \M + -X\ 2 M' + \X\ 2 M'_ 



< 0. 
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This achieves the proof. 



□ 

Next, we wish to study F(X) when A approaches the segment [M_, M + ], assuming that 
M G Cg. The idea is to use a well known result from the theory of Cauchy integrals (see 
[8], |9], and [10]). Let us first give some technical definitions. 

DEFINITION 3.4. Defining the complex open halfspaces as C ± = {A G C / ±lm A > 0} 

we say that h± : C ± — ► C is locally Holder continuous in C ± with exponent 7 G ]0, 1[ if 

V J R>0, 3C R >0 I V(Ai,A 2 )G^, M*i)-M*2)| <C r |Ai-A 2 | 7 (3.25) 

where B% = {A G C± / |A| < R}. If the constant C R can be taken independent of R we say 
that h± is Holder continuous in C*. 

In the same way we say that h± : \ {M_,M + } — ► C is locally Holder continuous in 
C± \ {M_,M + } with exponent 7 G ]0, 1[ if 

VR,5>0, 3C RiS >0 I V(Ai,A 2 ) e B% s , |/i±(Ai) - /i±(A 2 )| < C R>5 |A X - A 2 p 

(3.26) 

w/zere S± 5 = {A G / |A - M±| > 5}. 

We quote below a classical result in Cauchy integrals that we use. 

PROPOSITION 3.5. (Plemelj-Privalov Theorem) Let g{z) : [M_,M+] ^ C be Holder 
continuous with exponent 7 G ]0, 1[ in [M_,M + ]. Define the function 

5(A) =/ g(z) (z - \y l dz, XeC\[M_,M + }. 
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Then, for any A G C\]M_, M 4 



Ma 



lim 5"(A ± is) — S±(X) :— P.V. I g{z) (z - A)" 1 dz ±ing(\), (3.27) 
£ i° Jm- 

where P.V. stands for the Cauchy principal value of the integral. The function 

S±(X) := S(X), !/A6C±\[M_,M + ] and S±(X), ifXe}M_,M + [ } (3.28) 



is analytic in C ± and locally Holder continuous with exponent 7 in C ± \ {M_, M + }. More- 
over, if g(M + ) = g(M-) = , S±(A) can be continuously extended to mto a function 
which is locally Holder continuous in C ± . Conversely, if g(M±) 7^ ; <S±(A) blows up as a 
logarithm when A tends to M±. 

In what follows, we assume that M G Cf' 7 . For the sequel, it will be useful to introduce 
fi(z) the inverse function to M(y), 



fi(z)=y, forz = M(y), y6 [-1,1]. 
Changing the variable of integration to z = M(y), we have that, 

fM+ 

F(X) = / fjf(z) (z - X)~ 2 dz, AgC\[M_,M 4 

JM- 



In order to apply Proposition 3.5, we integrate by parts in (|3.30|) and obtain that, 

Ma 



F(X) = -//(z) (z - A) 



-1 



+ / fj,"(z) (z - X)' 1 dz. 

M- Jm_ 



(3.29) 



(3.30) 



(3.31) 



Applying Proposition 3.5 with g = /i", we get 



lim ei0 F(X ± is) = F ± (X), V A G]M_, M+[, 



F±(A) :=-/i'(z) (s-A) 



-1 



M 4 



Af 4 



+ P.V. / //"(z) (z - A) -1 dz ±infi"(\), 

M- Jm_ 



(3.32) 



and we can define as in Proposition 3.5, formula (3.28), F±(X) as a function that is analytic 
in C ± and locally Holder continuous with exponent 7 in C ± \ {M_,M + }. Note that for 
A G ]M-,M+[,Im F±(A) = ±7r//"(A) ^ 0, in particular F±(A) -2^0, which allows us to 
state the following result. 
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LEMMA 3.6. Assume that M G C s 2 ' 7 . Then, for any A G ]M_, M 4 



lim N(X ± is) = (2 - F±(A)) \ (3.33) 



Moreover, the function 



N(X), A G = {z G C/ ± Imz > 0}, 

7V ± (A) := <( (3.34) 
(F±(A)-2) \ AG]M_,M + [, 

«5 extended by continuity to C ± , wift iV±(M) = /or M = M + or M_ ; as a function 
analytic in C ± and Holder continuous in C ± w^A exponent 7. Moreover, 



VAG[M_,M+], JV_(A) = iV+(A), i e. JV±(A) = fte iV(A) ± iTm iV(A), (3.35) 
where Tie N is the common value of Re N + and Tie iV_ and Xm N := Im N + = —Xm iV_. 
REMARK 3.7. One can easily obtain an explicit expression of Tie N and Xm N from 



(3.32 ) and (3.28). For instance, one has 



XmN{\) = 7i\N ± (\)\ 2 AG[M_,M+]. (3.36) 

Such expressions are, in particular, useful for numerical computations (see 



Proof: The continuous extension by N±{M±) = is valid by (3.31). The local Holder 
continuity of N± in \ {M_,M + } is inherited from the same property for F±. The 
Holder continuity up to M± remains to be clarified. We consider the case of M + (M_ is 
treated analogously) and introduce the function 

J+(A) = (M+ - A) (2 - F±(A)), such that iV±(A) = (M+ - A) J + (A) -1 . 

We can conclude if we show that J + is Holder continuous in a neighborhood of M + and 
does not vanish at M = M + . According to ( 3.31[ ), 
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that we rewrite as 



J+(A) = //(M 4 



iJI + -X) (2-^4-/(M + )lo^ 



- (M+-X) 



M_ - A 
M + /i»( z )-/(Jlf + ) 



M+ - A 
M_ - A 



(3.37) 



Af_ 



z- A 



dz. 



Consequently, since 



ihi logo; is Holder continuous with any exponent in (0, 1), 



by Proposition 3.5 applied with g{z) = /jl"(z) — /i"(M + ) (that satisfies g(M + ) = 0), 



A 



M- 



Z-X 



dz is Holder continuous with exponent 7, 



we deduce from (3.37) that J+ is Holder continuous with exponent 7 in a neighborhood of 



M+. Finally, J+(M + ) = fi'(M + ) ^ 0, since M e C 2 ' 7 . 



□ 



We shall need an analogous result for N'(X), the derivative of N(X). For this, we assume 
that M e C s 3 - 7 and remark that for A G C \ [M_, M+] 



F'(A) 



Ai 4 



Ai_ 



dz 



/j, {z) — [(z - A) 2 ] <iz 



(z - A)- 2 



M. 



+ 



+ 



M- JM 
M+ 



fi"{z) {z-X)- 2 dz. 



M- 



p!\z) (z - A)- 1 



AU 



I n'"(z) (z - A)" 1 dz. 

M_ J M _ 



(3.38) 



Applying Proposition 3.5 as for F(X), we obtain 



f lim ei0 F'(\±ie) = F4(A), VAe]M_,M+[, 



F4(A) := -//(^-A)- 2 -^(^-A)- 1 



M+ 



Ai_ 



(3.39) 



Ai 1 



P.V. / n"'(z)(z-X)- 1 dz ±V(A), 
'm_ 
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and we can define again as in Proposition 3.5, formula (3.28), F±(\) as a function that is 
analytic in C ± and locally Holder continuous with exponent 7 in C ± \ {M_,M + }. Next 



we state for N' the equivalent result of Lemma [373] for N. 
LEMMA 3.8. Assume that M G C s 3 ' 7 . Then, for any X G ]M_, M 4 

\unN'(X±ie) = -F^X) (2-F ± (A))~ 2 . 



(3.40) 



A G = {z G C / ± Im z > 0}, 



Moreover, the function 

r ;v'(A), 

iVi(A) := 

( -F4(A) (2-F ± (A))~ 2 , AG]M_,M + [, 
is extended by continuity to 

N'AM) = - ?L±—L forM = M + orM_. 
±v ; jj,'(M) 2 J + 



(3.41) 



(3.42) 



TTien, N±(X) is analytic in C ± and Holder continuous in with exponent 7. Moreover, 



VAe[M_,M + ], Ni(A) = jV|(A), ». e. iV"4(A) = fte JV'(A) ± iXm JV'(A), (3.43) 
where Tie N' is the common value of Tie N' + and Tie N'_ andlm N' = Im N' + = -Im N'_. 



Proof: It is very similar to the proof of Lemma 3.3[ From formula (3.38), using the same 



trick used for proving the Holder continuity of J + (X) (see the proof of Lemma 3.3), we 
deduce that: 

F4(A) = {M+ - X)- 2 [ - //'(M+) + h{\)] , h G C ' 7 , h(0) = 0. 
Therefore, since iV±(A) 2 = (M+ - A) 2 J+(A)- 2 , 

N' ± (X) = -F4(A) iV ± (A) 2 = [ - //(M+) + h(X)] J+(A)- 2 . 
It is then easy to conclude. 



□ 
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REMARK 3.9. Again we can easily obtain an explicit expression of Tie N' and Tm N' 
from p39p and ( pT34p (see 



With Lemmata 3J3 and 3.8 , we now have all the needed information for studying the 
integrals Ii(kt,y),£ = 0, 1. We first obtain a new expression for Ii tC (kt,y): 

LEMMA 3.10. Assume that M e , then 

[iV + (A)-AL(A)] 



If moreover, M e Cf ,7 ; then 



AU 



P.V. I < 

<M. 



-ikXt 



dX 



(A - M{y)) 
m[N + (M{y)) + JV_(Af(y))] e - ikM{y)t . 



(3.44) 



h, c {kt,y) := -iktI 0iC (kt,y)+P.V. 



AU 



M- 



(A - M(y)) 



(3.45) 



- m [ N' + (M(y)) + N'_ (M(y)) ] e" ikM ^. 



IoAkt,y) := 



± *7? 



Proof: Let us start from the formula (3.22) for I 0lC (kt,y) that we can rewrite: 

N ±M (X) ^ f M+ N kt (X±i6) 

(A - M(y)) Y' 

where N± M (X) := N±(X) e~ iktx . Since N(M±) = 0, 

N kt (X) 



M _ (X±i5-M(y)) 



dX 



lim 



6-* J lf (X-M(y)) 



dX = 0, 



so that 



s-^Jm- {X + i5-M(y)) ^oJ M _ (\-i8-M(y)) 



that we rewrite as 



I , c (kt,y) = \l&\kt,y)-lfc-\kt,y)] + \ I& + \kt,y) - I&-\kt,y) 



where 



<5^0 



M + N± M (X±i5)-N ±M (X) 



lgt\kt,y) = lim / 



(A±i«J-M(3/)) 
W±,fct(A) 



w (X±i6-M(y)) 
21 



The Holder continuity of N±^ kt (\) (see Lemma 3.3) implies ^^(kt.y) = 0. Then, it is 
sufficient to apply again Proposition 3.5 with g = —N(X) e~ lkt to get the result. 



For Ii tC (kt, y), we first integrate by parts in ( |3.22[ ) to obtain 

h, c (kt,y) := 



dX. 



, Ss (X-M(y)) 

We then proceed as above using N± kt (X) = (N±(X) — iktN±(X)) e~ lktx and Lemma 
The details are left to the reader. 



3.8 



□ 



We are now in position to derive our estimates for Ii(kt,y),£ = 0, 1. 

LEMMA 3.11. Assume that M £ C e s +1 ,£ = 1,2. Then, there exists a constant 2| > 0, 
depending only on £ and M, such that one has the uniform estimates: 



WkeR, Vt>0, \\I e (kt,-)\\ L oo i{ _ 1;1]) <I^ (l + £\ktf), £ = 0,1. 



(3.46) 



Proof: Remarking that, thanks to (3.35), (N+ + iV_)(A) = 2TZe N(X) and setting 



AiV(A) := (N+ - AL)(A) = 2i(lm N(X) = 2iir \N ± {X)\ V'(A), 



(3.47) 



by (3.35) and (3.36). From (3.44) we deduce that, 



V(Art,y) <\P-V. / 



« tfcAt ,. AiV ^A, dA|+27r7eeiV(M(y))|. 



(A - M(y)) 



(3.48) 



By Lemma 3.3 we know that AiV(A) is Holder continuous with exponent 7. Since AiV(A) = 
0, for A < M_ and for A > M+, we can write, denoting L = M + — M_ (so that M(y) — L < 
M_ and M(y) + L > M+), 



P.V. 



M+ 



-ik\t 



M- 



AN(X) 
!A-M(y)) 



dA = P.V. 



M(y)+L 
M(y)-L 



(A - MO,)) 



e -ikM(y)t p y 



L 
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V 



(3.49) 



Next, we write 



P.V. 



V 



-ikvt 



AN(v + M(y)) - AN(M(y)) 



+ AN(M(y)) [P.V. 



v 

-ikvt 



dv 



v 



Finally, we remark that 

/L —ikvt 
dv -- 
-L " 

is bounded (in modulus) since 

30 sin(0 



sm(kvt) , 
-1 I — dv 



* LU sin(0 

Lkt S 



Finally, if we set 



K : = sup 

A>0 



sin(Q 



lim 



sin(0 



< +00. 



and I AiV| T := sup 



we deduce from (|3.49D and (|3.50|) that 

-L 



P.V. I ^ < 



-ikvt 



AN(v + M(y)) 



dv 



v 



< 



L oo meaning 
L dv 

L 



\v\ l -~i 



\AN{x) - AN{y)\ 
\x-y\i 

■ \\l°°(-M-,M+))' 

|AiV| 7 + K \\AN\\ L . 



which, using (3.48), gives: 



|io, c (*t,y)| </, 



O.c 



2 V^ ^- \ANL + K \\AN\\ Laa + 2?rptiV L <». 
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Finally, since (see (3.21, 3.22)) 



Mkt,y)\ <2nJ2 | A R "_^| + \W(kt,y)\ 



we obtain the inequality (3.46) for £ = with 



/* = /* c + 27r^|ResiV(A ± ; 



'A±-MQ 



dv 



(3.50) 



(3.51) 



(3.52) 



(3.53) 



The inequality (3.46) for £ = 1 is derived analogously using (3.45) and Lemma 3.8 This 



time, instead of AiV, we have to work with AiV' := N', — N'_ which fortunately satisfies 
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AN'(X) = for A G (-00, M_] U [M + ,oo), thanks to (|3.42[). Note that the additional 



factor \kt\ comes from the first term in the right-hand side of (3.45). 



□ 



We now go back to the estimates of a and k 1 . Using the fact that |a(w)| < ||w||z,2 \\v\\ L 2 



(Cauchy-Schwartz) and \\u\\l 2 < V% \\u\\l^, we deduce from (3.13) and Lemma 
for n — 0, 1, • • • , 

\k\ 1+n |a (M)| < ^ {l*o +I{ \\M\\ L ~{l + \kt\)) \k\ 1+n ||«°(MIUg- 



3.11 



that 



(3.54) 



In the same way, from (3.14) and again Lemma 3.11, we obtain 



/2 

A;| 1+rt |ai(M)| < (l + M) \k\ n llw^OlUg, n = 0,l,-- - . 



(3.55) 



Thus, by Plancherel's theorem, we see that, C denoting a constant depending on ||M|| LOO , I£ 
and II (and thus on M only), 



d 1+n a(u) f d 1+n u° ii 9 ^ 1 !! 

dx l+n V'Vm ^ 6 { II 5x l+„ him + \\-g^r\\Ll(Ll) 



+ t 



d 2+n u° , 



Q x 2+r 



Ql+n„.l 



(3.56) 



Q x l+r 



for n = 0, 1, • • • . 



From (3.1), or equivalently (3.3), one easily deduces that for each y G] — 1, 1[ and t > 0: 



/ 5 
y, t) = u°{x - tM(y), y) + t \il\x - M(y)t, y) + M{y)—u°{x - M(y)t, y) 

+ jl dzi J Z1 dz 2 ^{x- tM(y) + z 2 M(y),z 2 ), 



(3.57) 



for all t > 0,(x,y) 6lx [—1,1]- Taking the L 2 -norm in y and using (3.56), we easily 
obtain (with another constant C depending only on M) 
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d n u(x,y,t) 
dx n 

\M\\ L o 



n n ,0 



\Ll(Ll) 



< 



O n v! 
dx n 



d n u (-, 



d 1+n u° 
dx 1+n ( 



+ Ct 2 ( 



5 l+n M l 



d 2+n u° 
Qx 2+n " L y( L ^ ~^ 1 1 dx 1+n 



(3.58) 



<9x 3+ 



n 



0,1, 



9m 



Concerning — we observe that taking the time derivative amounts to the following: in 

at 

the first equation in (3.22) to multiplication by —ik\±, and in equations (3.44, 3.45) to 
multiplication by —ikX,X e [M_,M + ] and by —ikM{y). We easily see that arguing as 



, da. 



above one obtains estimates for \—-(k,t)\ and for |— — (k,t)\ which are similar to those 



Of 



Of 



for \& (k,t)\ and \a.i(k, t)\ up to an extra power of \k\. As a consequence, it is easy to 
obtain the following estimate for the mean value of a(u): 



t 



3l{U) 

dtdx n 

\\d 2 + n u° 



[X, 



Li 



< c 



^ t "m u || _ || d^-^u 1 || 

Q x 2+n ||L2(i|) -T || Q x l + n II -Lg (ig ) 



d 1+n u° 



dx 



d n u l 

l+n H^yC^i) II Q x n \\ L y( L V) 



n 



o,i,... . 



(3.59) 



Furthermore, taking the time derivative in both sides of (3.57), we obtain that, 



U 



dtdx 
Ct 



(x,y,t)\\ L2{Ll) < \\(1 + M(y)) 



d n+1 u° 



d n u l 



d 2+n u° || J_ || l9"+ 1 M 1 || I 

a^r||L2(L2) + \\-q^+t\\lI{LI) ) + 



dx n+T\\ L K L l) + I' -Q-^\\LliLl) + 

Cr \ || - - Wmhi) + 



<9x 3+r 



<9x 2H 



\Ll{Ll) 



| <9 4 +"-u° || , || d 3 + n u x || i 



0,1, 



(3.60) 



In the proof of estimates (3.58, 3.59, 3.60) we were not bothered by factor 1/k in the 



right-hand side of equation (3.14) because we only needed to estimate 9 + dx 2 U ^ what implies 
that we have to multiply a(w) by (ik) 1+n , what removes the singularity at k = 0. However, 
as a(w) is a relevant physical variable it is important to estimate it directly. To do this we 
proceed below in a slightly different way. We define, 

1 



I(k,t,y) :-- 



2nk 
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h(kt,y). 



(3.61) 



Then, equation (3.14) is replaced by, 



an(k,t) = -2a(l(k,t,-)u 1 (k,-)). 



We regularize the singularity at k = in / as follows. By Cauchy's theorem, 



N(X) 



and then, 



It follows that, 



Xl (A - M{y)f 
N(X) 



dX 



N(X) 



R+iAj 



(X-M(y)y 



dX = lim 



r (A-M(y))> 

N(X) 



dX, for R > Xj, 



dX = 0. 



dX, A/ > 0. 



2tt J R+lXi k (A - M{y)f 
This representation of I is useful because ( e ~* fcAt — is not singular at k = 0. 

Deforming the contour of integration as above, we prove that, 



where 



and 



I p := -z- 



/ = / p + J c , 



1 ResMA_) e -ik\+t_ l R es iv(A^ 
z- 



k {X_-M{y)Y k (X + -M(y)Y 



M+ 

I r = P.V. I dX 

'M- 



t e -*At|jV(A + iO)|V(A) 



A - M (y) 



+ i {e -i k xt_ 1) c sNl{x + i0) i 



(A-Mfo)) 



dA 



-te- ikM ^lR.N(M(y) +i0) - | (e~ ifcM ^ - l)3t.N'(N(y) + iO). 



Note that as, 



we have that, 



/A' 



A< 



/A- 



<• —ikXt -p 



< At. 
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As in the proof of Lemma 3.11 we prove that 



VfceR, Vt>0, ||/(fc,t,-)|U~([-i,i]) < I* (! + *)• 



(3.69) 



Arguing as above and using (3.9, 3.13, 3.14) Lemma 3.11 and (3.61 3.69) we prove that 



|a(«)|U<C 



kt u +t 



du c 



dx 



L 2 y(Ll) 



+ (l + t) \\u l \ 



L 2 y(Ll) 



(3.70) 



Let us state the result above theorem: 



THEOREM 3.12. Let us assume that M G Cf' 7 . Then, Cauchy problem (1.1) is weakly 



well-posed in the sense that for any (vP^u 1 ) G L 2 y (Hl +n ) x L 2 y (Hl +n ),n = 0, 1. • • • , it 
admits a unique solution satisfying 



(3.71) 



and the estimates (3.58, 3.59, 3.60 j and (3.70) hold. 



The estimates (3.58, 3.59, 3.60, 3.70) deserve some comments. 



Concerning the mean value of the solution a(w), the estimate (3.70) predicts: 



A loss of regularity in x (by one order of regularity) between the initial data 
and the solution at time t > 0. 

A corresponding polynomial growth in time (t) of the solution. This phe- 
nomenon is similar to the one observed with weakly hyperbolic systems ( seefTT]). 
Such a phenomenon does not occur, of course, in the case of a uniform reference 
flow but we think that these results are sharp in our case. 



Concerning the solution itself the estimate (3.58) announces a possible loss of regu 



larity in x by three orders, together with a corresponding polynomial growth in time 
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as t 3 . We do not know if this estimate is optimal: the way it has been obtained is by 
solving the square transport equation but ignoring that the right hand side depends 
on a.(u), i. e., precisely on the solution we are looking for. It might be that the 
additional lost of two derivatives is an artefact introduced by our technique. 

3.3.2 For piecewise linear profiles 

We consider now profiles M(y) that are continuous, piece-wise linear and strictly increasing 
or strictly decreasing. That is to say, such that there are finite sequences, —1 = x_ N < 
X-n+i < ■ ■ ■ < x < xi < ■ ■ ■ < x N = I and M, < M i+1 , or M { > M i+1 , i e {—N, -N + 
l,-- - ,N-1} with, 

V-N<i<N-l, M{y) = a i y + i3 i ,y e [x h x i+1 ], 

M i+1 -Mj 

By eventually redefining the partition we can always assume that a i+ i ^ oij 



a i = s-TI-s-S a i x i + A = Mi, I = i, i + 1. 



(3.72) 



By explicit computation we prove that, 

N-l 



M i+ i - Mi 

i=-N 1+1 1 



(3.73) 



Mi - X M i+1 - A 

We see that F(X) extends to a meromorphic function on C and that the equation F(X) — 2 
is a polynomial equation of order at most 2N (there can be cancelations). 

We make now the assumption, 



(Hsl) N(X) has no complex (i.e. non real) poles, and in [M_,M+] it has L poles 
that are all simple. 

Note that since N(X) has the poles A_ < M_ and A + > M + necessarily L < 2N — 2. In 
[3] it is proven that TCsl holds if M(y) is odd increasing and convex or odd decreasing and 
concave. In both cases L = 2N — 2. Furthermore, [5] proves that Hsl holds if M(y) is 
even, increasing or decreasing, and either convex or concave. In all these cases L = N — 1. 
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Suppose that the poles of N(\) in [M_, M + ] are located at the points \j with residues 



rj,j = 1,2, 



L. 



Equations (3.9, 3.13, 3.14, 3.22, 3.61, 3.62, 3.64, 3.65) remain true in this case. The 



calculation of Ii c ,l = 1,2 in (3.22) and of I c in (3.64) is much simpler now. By the residues 



theorem it is given by the sum of the residues of the integrand at the poles Xj,j = 1, 2, ■ ■ • , L 
and M(y). We have that. 



I 0)C = -27U 



X, - M(y) 



E 



+ e~ ikM ^ t N{M{y)) 



Ji )C = -2m 



E 



g ikXjt 



- (\j-M{y)y 



ikt e~ %kM{v)t N(M(y)) + e ~ ikM{y)t N'(M(y)) 



(3.74) 



(3.75) 



E 

3=1 



-ik (Xj - M(y)¥ 



te~ ikM{y)t N{M{y))- 



(3.76) 



{„ — ikM(y)t i \ , , . . „ 

{ - ^±N'{M{y)). 



-ik 



We prove the following estimates reordering the terms in Ij )C , I = 1,2, J Cj , and developing 
e ikt(\j-M(y)) j n a Baylor expansion when ikt(\j — M(y)),j = 1, 2, • • • , L is small. 



IVI <C(l + \k\t), 



0_ 

Ft 



< C\k\{l + \k\t), 



(3.77) 
(3.78) 



\h, c \ <C(1 + (H) 2 ), 



d_ 

Of 



< C\k\{l + {ktf 



\I C \ < C(t+\k\t 2 ). 



(3.79) 
(3.80) 
(3.81) 



Inverting the Fourier transform in (3.9) and using (3.13, 3.14), the first equation in (3.22) 
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and (3.77 3.80) we prove that, 



,d 1+n a(V) , „/ „d 1+n u° 



.. _/ ..c^u" n^V., n d 2+n u° 

dxl+n c^tjikg < c { \\-Q^\\mLi) + w^uiiLD - t 



Q x l+n u^v 

d 3 + n M° II I || a 2 +"« 1 n 1 \ 



71 = 0,1,-- 



d n u(x, y, t) 



i(x,y,t) d n u° ( d n u x {-,y) 

-fan ^ W-Q^rvMli(^)+ t { II dx n + ll M II 



(3.82) 



d 1+ "u° / „,\|| II d 2+n u° || i n a^M 1 !! \ j_ 3 I, d 3+n u° 

-Q^r[-,y)\\Ll(Li) J + ^ [t ^ Hg^+rrll^^l) + \\-Q^+n\\Ll(Ll)) + * II H 



+t 4 



I 8 4 +"m° [I _| || d 3 + n y} || 

I fe 4+u "t" II Q x 3+ n \\Ll(Ll) 



71=0,1, 



(3.83) 



d n v} d 2+n u° 

^rll^(^) + * II Q x 2+n II 



n = 0,1, 



(3.84) 



-(x,y,t)\\ Ll(Li) < \\(l + M{y))—— || LS(L|) + ||__ \\ L 2 (Li) + 



dtdx r 



dx r - 

d 3+n u° 

Q x 3+n H- L y( L l)~'~ 



- ( \\ d 2+n u° || , || g"+V || \ , + 2 ( 

\\-q^+^\\lI{lI)+ II -Q^+r || Ll{Ll)J + t ^ 

a 2+n -u 1 || A i +3 f ||(9 3 + n « 1 || i ||d 4 +'V|| \ i 

-Q^\\L%{Ll) J + t (J| -g^r\\Ll(Ll) + ||^3+^||L2(L2)j + 



+i( || d 5 + n u° || | n^+^M 1 !! 

1 II a^5+7T||L2(L2) + ||^4+^||L2 (L 2) 



71= 0,1, 



(3.85) 



and using also (3.81 ) 



|a(tt)||„<G 



\ a \\Ll(Ll)^ 1 



du° 



dx 



L 2 (L 2 ) 



+ t IKIlL2(L2) + 



d 2 u° 




du 1 




dx 2 


+ 


dx 


L 2 y(L 2 )). 



(3.86) 



Let us state the result above theorem: 
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THEOREM 3.13. Let us assume that M(y) is strictly increasing or strictly decreasing 



and that it satisfies assumption Hsl- Then, Cauchy problem (1.1 j is weakly well-posed 



in the sense that for any (u°, u 1 ) G Ly(H^ +n ) x L?(Hl +n ),n — 0, 1. • • • , it admits a unique 
solution satisfying 

u e C\R + ; Ll{H n x )) n C\R + ; L^" 1 )), (3.87) 



and the estimates (3.82 -3.86 j hold. 



Note that in comparison with the case of smooth profiles we loose an extra order in 
regularity and we have an extra power of increase in time. We believe that this not just a 
consequence of our method, and that it is actually due to the discontinuity on the derivative 
of M. 

4 A quasi-explicit representation of the solution and its 
physical interpretation 

4.1 For a class of smooth profiles 



We consider that we are in the conditions of Subsection 3.3.1 Applying the inverse Fourier 



transform to (3.9), and using (3.13, 3.14) the first equation in (3.22), (3.44, 3.45) and (3.61 



3.62, 3.64 3.67) and using the fact that multiplication by e tks in Fourier space amounts to 



translation by — s in the physical x space, we obtain the following representation of &(u) 



&(u)(x, t) = 2L p (u)(x, t) + a c (w)(x, t), 



(4.1) 



where, 

a p (-u)(x, t) :- 

ResiV(A + ) 
ResiV(A_) 



2« { Io +t u\x-s,y)ds + u\x - s,y) ds- 



(A+-M( 3/ ))2 (\+-M(y)) 



X- 



u°(x - X+t,y)- 



u°(x — A__t, y) 



(4.2) 
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is the contribution of the poles, A±, of N(X) (or equivalently of the point spectrum of the 
operator A), and 

a c (u)(x,t) = -2a [p.V. /£+ d\ [ffi^V(A) (tu\x - Xt,y) + tM(y)£u°(x - \t,y)+ 

u°(x - At, y)) + { J* u\x - s, y) ds - M(y)u°(x - At, y) }" - 

.KeN(M(y) + iO) {tu l (x - M(y)t,y)+ tM(y)^u°(x - M(y)t,y) + u°(x - M(y)t,y)} - 
Tie N'(M(y) + iO) f ff (y)t v}{x - s, y) ds - M(y)u°(x - M(y)t, y) 



(4.3) 

is the contribution of the branch cut of N(X) (or equivalently of the continuous spectrum 
of A). 



Equations (3.57) and (4.1, 4.2, 4.3) give a quasi-explicit representation of the solution, 



u(x,y,t) to problem (1.1). 



Our quasi-explicit representation gives a nice physical interpretation for the propagation 
of acoustic waves in the fluid, that we present below. Suppose for simplicity that u 1 = 0. 



Equation (4.1) for the average of the solution can be written as follows, 



a(x, t) = a p (x, t) + P.V. / d\a.^(x,t, A) + / dX a^\x, A,t) 



-M 4 

' M- JM- 

where a p (x,t) is a solution of the generalized square transport equation, 



d B\ (d , 

dt +X+ dx) U + A -^ )ap(M) = ' 



(4.4) 



(4.5) 



ag\u)(x, A,t) := -2a 



Xm { Z { M^ M{y)u»{x-\t,y) 



is a solution of the square transport equation, 



| + aA] tf)(„)(x,A,t) = 0, 



(4.6) 



(4.7) 
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and, 

a( 2 )(x,A,t) :=-//(A) 



-Tie N(X + zO) jtA^u°(ac - At, //(A)) + - At, //(A)) J + 
fte JV'(A + zO)Am°(:t - At, /i(A))] , 
is also a solution of the square transport equation, 



(4- 



£ + a£) ^(.,a,*) = o. 



(4.9) 



In the case where w 1 7^ we have a similar representation for t| [a(w)j (x,t). 



4.2 For piecewise linear profiles 



We consider that we are in the conditions of Subsection 3.3.2 Applying the inverse Fourier 



transform to (3.9), and using (3.13, 3.14) the first equation in (3.22), (3.61 , 3.62, 3.64 3.67) 



and (3.74 3.76), and using the fact that multiplication by e lks in Fourier space amounts 



to translation by — s in the physical x space, we obtain the following representation of a(u) 



a.(u)(x, t) = 2L p (u)(x, t) + a. c (u)(x, t), 



(4.10) 



where, 



2Lp(u)(x, t) := 2a 



ResiV(A+) 
(A+ - M{y)f J 



x+t 



m 1 (x — s, y) ds+ 



ResiV(A_) 
(A_-M(y))2 

u°(x - X+t,y) - ResiV(A_) 



/•X-t 

/ u x (x — s, y) ds — Res iV(A- 
Jo 



(X--M(y))2 (X—M{y)) 



(A + - M{y)Y (A+ - M(y)) 
u°(x - A_t,y)J , 



(4.11) 



is the same as in (4.2), i.e., it is the contribution of the poles A± of N(X) (or equivalently 
of the point spectrum of the operator A), and a. c (u)(x,t) is the contribution of the poles 
of iV(A) in [M_,M + ] (note that [M_,M + ] is the continuous spectrum of the operator A), 
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and it is given by, 

L 

0(s, t) = 2a{M(y) [ - £ — _^ (y))2 u°(* - A,f, y) + tN(M(y)) 



slAu) 



£u°(x - M(y)t, y) - N'(M(y))u°(x - M(y)t, y) + £^ =1 - A,-*, y) 



(4.12) 



+iv(M(y)) M o (x-M( z/ )t,y) + + ^ =1 jx^kvw ti jtul ( x - s >y) ds 

+tN(M(y)) u\x - M{y)t, y) + N'(M(y)) J Q M(y)t u\x - s, y) ds]. 
As for smooth profiles we give a physical interpretation of the solution. For simplicity 



we consider only the case u l = 0. Equation (4.10) can be written as follows. 



a(x,t) = sl p (x, t) + / dXa. c (x,t, A) (4-13) 

JM- 

where, as for smooth profiles, a p (x,t) is a solution of the generalized square transport 
equation, 

and 



a c (x, t, A) := 2//(A) { [ - A ( A .^ A )2 - A i*> MA)) 

u°(x - Xd, u(X))] + tXN(X)^-u°(x - Xt, u(X)) ( 415 ) 
j — X J ox 



A, - A 

-XN'(X)u°(x - Xt, //(A)) + iV(A) m°(x - Xt, y) 
is a solution of the generalized transport equation, 



n(| + >4)(| + A|)>(^H0. ,4,6) 
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